Dynamic model of spherical perturbations in the 
Friedman universe. III. Automodel solutions. 
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Abstract 

A class of exact spherically symmetric perturbations of retarding auto- 
model solutions linearized around Friedman background of Einstein equa- 
tions for an ideal fluid with an arbitrary barotrope value is obtained and 
^ I investigated. 

00 ; 1 Introduction 

Q ' In the previous works of the authors [1,2] a class of exact retarding solutions 

for hnear spherical perturbations of the Friedman universe with an ultrarela- 
5^ ' tivistic equation of state of the ideal fluid filling it, corresponding to the central 

bi), singular source presence and having the sight of polinoms by radial variable, 

was obtained. At that it was noted, by zero boarding conditions at the sound 
horizon for the class metrics perturbations the energy dense perturbations 
^ ' have the first genus break at the sound horizon. In this paper we shall study 

OO , the retarding solutions in detail by extending the investigations range to the 

equations of state of the fluid with an arbitrary barotrope coefficient k. 

So we shall investigate the retarding solutions of the evolutionary equation 
for the spherical perturbations [1] (84): 



X 



•■ 2 • 6 1 -t- k) * „ 

T] (1 + 3k)2 ^ ^ 



with boarding conditions at the sound horizon corresponding to the zero values 
of the metrics perturbations and its first radial derivatives 

; So : 7" = V^V- (2) 

'^ir,v)\rv^, = ^^{v); *V,^)l.^, = o, (3) 

where /i(?7) is the central source singular mass, so that the metrics component 
perturbation 944 is equal to 

Sg,,^a'iv)S., s.^2^ = 2^^^!:i^!l^J^, (4) 

ar ar 

moreover the function ^(r, rj) corresponds to the nonsingular part of the poten- 
tial: 

*(0,?7) = 0. (5) 
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Temporal mass evolution is described by the equation 
2. 6(1 + «) ^1 

?/ (1 + dKj^ 



which has the solution 



3(1 + k) 



^ = ^+7?i+3k +/i_7? 1+3K , (1 + k)^0; (7) 
(for detail see the previous works [1,2]). 

The potential function ^(r, 77) and the scalar ii{ri) completely determine the 
energy density perturbations and the fluid velocity 

' ^3^$-*"V (8) 



So Anra^eo \ a 

1 9 $ 

Anra^Eo dr r ' 

where £0(^7) is a non-perturbed energy density of the Friedman universe 

6(1+k) 2 6k 

Ea ^ rj i+3k ; a^T^i+SK; equ^ ^ rj i+3k . (10) 



2 Automodel solutions 



2.1 A general automodel solution 

So we shall look for solutions of the evolutionary equations for the perturbations 
with zero boarding conditions at the zero sound horizon Supposing that in 



where 



(11) 



(12) 



we come to the automodel solutions class and get a common differential equation 
for the function Ga{z)'. 



{1- z^)G\{z) + 2azG'4z) + 



6(1 + k) 



a{l + a) 



G(z)„=0. (13) 



.(l+3At)2 

The potential function 4'(r, 77) must be a combination of the private solutions 
^ir,v)=J2v''Go.iz). (14) 
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From ([3]) and ([7]) it follows that by a 7^ this combination can have two members 
only 

3(1 + k) 



^{r,r)) =G+{z)rii+3'^ +G-{z)t] i+3« , (15) 

at that in consequence of the boarding conditions ([3]) the functions G±{z) must 
satisfy the following boarding conditions 

G±(l)=/i±; G'±(1) = 0. (16) 

In particular from ^ and ([5]) it follows immediately that in the case a = 
the zero mass of the singular source corresponds to the pointed out class of 
solutions. 

At the parameter a arbitrary values the general solution of the linear differential 
equation is expressed by the Legendre functions, P^(z), and the adjoint 
Legendre functions Q'J^ 



G{z) = 



1 - z 



2(1+k) 2(1+k) 



(17) 



2.2 An auto model solution with the particle-like source 

The mentioned above is just for the formal general solution (|17|) at an arbitrary 
value of the parameter a. However, in our particular case (|15|) the values of the 
parameter a: 



3(1 + k) 



(18) 



1 + 3k ' 1 + 3k 
are the quadratic equation roots simultaneously 

^-o(lH-o) = 0. (19) 

Therefore in this case the equation (fT3| degenerates into a simpler one 

[l- z^)G\{z) + 2azG'a{z) = Q, (20) 

However, integrating it we obtain the following 

G'„(z) = Ci(l-z2)". (21) 



Comparing the second boarding condition (|T6l) with the expression ([21]) we 
see that for the fulfillment of the zero conditions for the first derivative of the 
potential at the zero sound horizon it is necessary that a > 0, that is The 



^ For example, see [7] and [8]. 
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particle-like source mass is growing by the time and is equal to zero at the time 
moment rj = 0; thus, =1 to provide a smooth sewing together of the solution 
with the Friedman one at the zero sound horizon it is necessary that 



Ai- = 0. 



(22) 



In the case of a > Oin consequence of the relation (|2T|) the second boarding 
condition is fulfilled by fulfilling the first one automatically, thus, the con- 
dition a > fulfillment provides a smooth sewing together of the solution in 
the class at the sound horizon. 

Now integrating the equation (j2ip formally and considering the conditions at 
the beginning of the coordinates ([5]), according to which 

G(0) = 0, (23) 

we find its formal solution within the whole interval of the values r — [0, -foo): 

i+3^+^a;dx, (z > 1), 



G{k,z) = Ci <^ 



2T 



1+3k 



where F(a,b,c,x) is a hyper-geometrical function (for example ,see [7]): 



r(/j)r(7 - /3) 7o 

3?(7) > n{l3) > 0; |arg(l - z)| < tt. 
At that the useful limitary relation is just [7]: 

Jm^F(a,/3,,,.).£|2)I^^^^ 5i(, - a - /?) > 0; (26) 

where r{x) is F-function. In Fig. 1 the solutions p4)) for values series of the 
barotrope coefficient are shown. 
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Fig. 1. Normalized function G{k,z) calculated with the formula 
(p4)) at the value ((27)) of the constant Ci and /i = 1. In the left part 
of the figure from bottom up k = 0; ^; ^; |; 1. 



In terms of the obtained solution ([24)) it is easy to see G{—z) — G{z), that is 
the function \['(r, 77) is an odd function of a radial variable really, it was shown 
in Ref. [2] . Considering the first boundary condition (IT6l) , we find the constant 
Ci: 




and thus, finally we obtain the automodel C^class solution corresponding to the 
zero boundary conditions at the sound horizon 



2 

^{r, 77, k) = ^+7] i+3« X 



2r(TW + i) r fl 2 3 r2 



^r(^ + l)V^^ V2' 1 + 3k'2' KT^^J' ^ (28) 
1, (r > \/Krj), 

Taking into consideration the relation according to 1^ we obtain the scalar 
function of the metrics perturbation 6iy: 
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5iy{r, 77) — — X 
r 



2r 



1 - 



\^ 1+3k 2 J 



L+.3K 



where x(x) is the Heavyside function: 
X{x) = {0,__a; 0;l,._a; > 0. 



x{V^V-r), (29) 



(30) 



In the private values series of the barotrope coefficient k the obtained solution 
is expressed in elementary functions k = 1/3 the ultra- relativistic fluid 



G(l/3,z) 
3 



1 3 

z z-^ 

3 



-u+ri z 



2^^^ 



1 - 



3 r 



1 / 



2V3?7 2 V%/3r/ 



(31) 



Being auto-model the solution ([5T|) coincides with the general solution in the 
form of the degree series, earlier obtained in the works [3,4,5,6], it confirms the 
correctness of the proved in [2] theorem about the retarding solution uniqueness 
in the case of the ultra- relativistic state equation once more, k = 1: fluid with 
extremely stiff state equation 



G(l, z) = —z\/l — + — arcsinz; 



TT \ 77 



7]^ 



2 arcsin — 



(32) 



Note, the first derivative conversion by a radial variable at the zero sound 
horizon in the obtained solutions at a > guaranties the relation (|2ip . 



2.3 A solution without a particle-like source 

Now supposing that in (fTTj) a = 0, we get the following equation instead of ([TF 



(1 - z^)G"{z) + 



6(1 -I-k) 
(1 + 3k)2 



G{z) 



0. 



(33) 



The general solution of this equation is a linear combination of the hyper- 
geometrical functions 
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G{k,z) = (1 - z^) X 



3k 5 + 9k 3 2\ ^ / 2 + 3k 3k — 1 1 2^ 



If the geometrical functions 

3k 5 + 9k 3 n\ , „ / 2 + 3k 3k — 1 1 , 



(34) 



, , z and F , , z 

V1 + 3k 2 + 6k'2' J Vl + 3« 2 + 6k'2' 

remained final at the sound horizon z = 1, the general solution ([55]) would auto- 
matically turn into zero at the sound horizon, then the zero boundary conditions 
at the sound horizon ([T51) would be in principle satisfied by the choice of con- 
stants in the general solution ([33|. However the pointed out hyper-geometrical 
functions have some peculiarities at the sound horizon. 

Really, as it is easy to see, the parameters a, /3, 7 of the hyper-geometrical func- 
tions of this linear combination satisfy the condition (See the hyper-geometrical 
function definition (26).) 

a + /3-7 = l, (35) 

At that for the first member of this combination the additional condition is 
fulfilled 

at the same time for the second one it is 

7.-/3.^+^. (37) 

In this case it is necessary to use the functional relation for the hyper-geometrical 
function [7] , which conformably to (j35p is given with the formula 

F{a, (3,a + (3 -l,z) = --^—-F{a - I, /3 - l,a + P ~~ 1)N (38) 
(1 — zj 



Using this relation in the formula (|34|) let us rewrite the general solution in the 
form 



G(k,z) ^(l + z) 



zCiF 
C2F 



1 



3(1 -hK) 1 ^2 



1 + 3k' 2(1 + 3k)' 2 

1 3(1 + k) 1 ; 

14- 3k' 2(1 + 3k)'~2''^ 



(39) 
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Calculating the hyper-geometrical functions values in the right part of Eq. ([5^ 
at the sound horizon we see these functions have peculiarities at the sound 
horizon. Therefore the auto-model solution satisfying the zero boundary condi- 
tions at the sound horizon in the case under investigation is the trivial solution 
Ci = C2 = only. Summarizing this subsection let us formulate the theorem 
Theorem. There are no retarding spherically- symmetric auto-model solutions 
of the equation ([T]) satisfying the zero boundary conditions (|16p at the sound 
horizon without the central particle-like source (fi = 0). 

The proved theorem is analogues to the theorem about the analytical solution 
of the Laplace equation in the spherical symmetry case. 



3 Research of the auto-model solutions 



3.1 Derivatives of the potential functions 

Let us come over to the analyses of the obtained auto-model solutions in the 
particle-like source case. For that let us use the expression (|S]) of the relative 
energy density of perturbation and the expression ^ of the radial medium 
velocity in the perturbation. At that we shall need the expressions for the first 
and second derivatives of the potential functions. Considering the definitions 
(fTT|) and (|12p and the relations ([20l) and (f27|) we shall get expressions for the 
first and second radial derivatives of the potential functions 'i'{r, r/) 



i-3k 2r 



1-I-3K 



kttF 



^"rr{r,v) = - 



8^+77 1+3k 



(1 - z2)l+3«^ 
f^+3\ 



(40) 



1-3k 



k{1 -\- 3k 



)^(t 



(l-z^)i+3«, (41) 



and for the temporal derivative of the function $(r, r/) also 
,G(fc,z)-/i+ 2/i+r 



$(r, 77) = 77H 



1+3k 



1 + 3k 



2 

L+3k 



1 



zil-z^) 



1+3k 



,(42) 



here the function G(k, z) is determined by the relation (|24l) with the constant 
Ci from (i27| . 

From the given expressions it follows: by n > —1/3 the first radial and 
temporal derivatives of the potential functions 4'(r, 77) and <f>(r, 77) turn into 
zero at the sound horizon 



— *(r,77) 



— $(r,77) 



= 0; (1 + 3k^0); 



(43) 
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0; (1 + 3k^0). 



(44) 



At K < 1/3 the second radial derivatives of the potential functions turn into 
zero at the sound horizon 



(92 



-■\/kr] 



= 0; (1-3k;5 0). 



(45) 



By K = 1/3 the second radial derivatives of the potential functions have a break 
of the first genus at the sound horizon. By k > 1/3 the second radial derivatives 
of the potential functions have a break of the second genus at the sound horizon. 



3.2 Evolution of the energy density distribution in the 
spherical perturbation 

Calculating the relative density of the perturbation energy by the formula ([8|) 
considering the relations (jlOp and (|40l) . (j4T|) we get finally 

5e 1 
£0 



Z777ry^(l + 3k) 



, G{k,z) - 11+ 
1 + 3k 



- 3- 



n V ( 7+9k \ 



3(1 + k) 
1+3k 



z{l - Z-^)1+3k + 



2/x+r 



7+9k 
2(1+3k) 



+3k 



1 



(l-z^) 



1-3k 
1+3k 



A(z)x(l-z), 



(46) 



where the reduced relative density of the perturbation energy A(z) is introduced. 
It can be strictly shown that A(z) ^ 0. 

From this expression it is seen that by the time the profile form of the 
perturbation energy density does not change relatively the dimensionless radial 
variable z = r j^fkr] and the relative density of the perturbation energy decreases 
in inverse proportion to the temporal variable 77 (Fig. 2). 

At that in the terms of the common radial variable r the energy density 
perturbation profile is deformed. For example in Fig. 3 the evolution of the 
relative energy density perturbation profile at the barotrope index k = 1/6 
is shown. =3 Further on from the formula ()46p it is immediately seen that 
at K < 1/3 the energy density perturbation at the sound horizon vanishes, at 
K = 1/3 it has a final jump at the sound horizon, and at k > 1/3 it has an 
infinite jump, it absolutely corresponds to the behavior of the second radial 
derivatives of the potential functions. 



9 



Fig. 3. It is the evolution of 
the reduced relative density of the 
energy A(z) calculated by the for- 
mula P5| as the function r at k = 
1/6. From the left to the right 77 = 
1;2;4;8;12;16. 



Fig. 2. It is the evolution of the 
reduced relative density of the per- 
turbation energy A(z) calculated by 
the formula (l45t as the function z. 
From bottom up in fine line k = 
1/6; 1/5; 1/4; 1/3; the heavy line cor- 
responds to K = 1/2. 

Let us find out the physical sense of the obtained solution. The energy perturba- 
tion corresponding to the nonsingular part of the potential function is described 
by the formula 

&E = ^Tic? / Ser^dr. 
Jo 

Hence considering (fTO|) we get 

SE = Atttj 1+3" / —r'^dr. 

Jo ^0 

Coming over to the dimensionless variable z in the integral by the formula 

r = \/7irjz, 

we bring it to the form 

2 /"i 

SE = Anfi+riWK K^/'^ A(z)dz -771(77), 
Jo 

where (see (7)): 



(47) 



777(77) ~ /i_(_77l+3K 



(48) 



is the mass of the central singular particle-like source. Thus, the full energy 
in the included in the nonsingular mode of the perturbation is negative and 
proportional to the mass of the particle-like source. 
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3.3 Evolution of the fluid radial velocity in the spherical 
perturbation 

Fulfilling the analogues calculations we get the expression for the radial velocity 
from © 



871772^3/2(1 + 3k) 



G{k,z)-^+ (1-z2)1+3k 



,, -n / 7+9k 
'^+-^ \^2(1+3k) 



1 



(l-z^) 



1-3k 
1+3k 



= 4t(z). 



(49) 



From this expression it is also seen that the radial velocity is negative and its 
profile remains constant within the scale z, and the absolute value of the velocity 
drops in inverse proportion to ij^. 




Fig. 4. Dependence of the profile of the perturbation radial 
velocity T(z) on the barotrop coefficient atK = 1/6. From bottom 
up K = 1/6; 1/4; 1/3; 1/2; 1 
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